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ABSTRACT 



O 



We construct a series of model galaxies in rotational equilibrium consisting of gas, stars, and a fixed dark matter (DM) halo and study 
how these equilibrium systems depend on the mass and form of the DM halo, gas temperature, non-thermal and rotation support 
^~», against gravity, and also on the redshift of galaxy formation. For every model galaxy we find the minimum gas mass M"™ required 

to achieve a state in which star formation (SF) is allowed according to contemporary SF criteria. The obtained M"™-Mdm relations 
are compared against the baryon-to-DM mass relation M b -M DM inferred from the ACDM theory and WMAP4 data. Our aim is to 
construct realistic initial models of dwarf galaxies (DGs), which take into account the gas self-gravity and can be used as a basis to 
study the dynamical and chemical evolution of DGs. Rotating equilibria are found by solving numerically the steady-state momentum 
equation for the gas component in the combined gravitational potential of gas, stars, and DM halo using a forward substitution 
procedure. We find that for a given M DM the value of M™ 11 depends crucially on the gas temperature T e , gas spin parameter a, degree 
of non-thermal support <x et y, and somewhat on the redshift for galaxy formation z g f. Depending on the actual values of T g , a, cr s g, and 
• Zgf, model galaxies may have M™" that are either greater or smaller than Mb. Galaxies with Mdm ~ 10 9 M Q are usually characterized 

by M™" S: Mb, implying that SF in such objects is a natural outcome as the required gas mass is consistent with what is available 
according to the ACDM theory. On the other hand, models with M DM Si 10 9 M Q are often characterized by M™" » M b , implying that 
they need much more gas than available to achieve a state in which SF is allowed. Our modeling suggests that a star-formation-allowed 
state is more difficult to achieve in DM halos with mass Si 10' M Q than in their upper-mass counterparts, because the required gas 
$H . mass often exceeds both Mb and Mdm- In the framework of the ACDM theory, this implies the existence of a critical DM halo mass 

ryj 1 below which the likelihood of star formation and hence the total stellar mass may drop substantially, in accordance with the stellar 

versus DM halo mass relations recently derived from the SDSS survey and Millennium Simulations. On the other hand, DGs that do 
not follow the ACDM trend are feasible and have recently been identified, which raises questions about the universality of the ACDM 
paradigm. 

Key words. Galaxies: dwarf - galaxies: structure - ISM: structure - Stars: formation - dark matter - Methods: numerical 
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' 1. Introduction viscous and diff erentially rotating) models of galaxies (e.g. 

iLindbloml |1992|) . It is therefore always realistic to assume that 
lO . The study of equilibrium states of self-gravitating, multi- the equilibrium configuration of a galaxy rotating about some 
O ■ component fluids is of considerable interest in astronomy be- axis is axisymmetric. It is however not always true that this fig- 
| cause they serve as basic models of many astrophysical objects ure of equilibrium is an ellipsoid. Rapidly rotating, compressible 
^ . (stars, protoplanetary disks, galaxies among others). While it gaS es characterized by a polytr opic equation of state quite natu- 
> ' is known (and quite obvious from simple symmetry consider- rally develop a flared structure (iBodenheimer & Ostrikerl Il973t 
ations) that isolated, non-rotating, self-gravitating fluids of finite |Tassoull[l978h. Flaring gas dis t ributions in some DGs have been 
extent must be spherically symmetric, it has been a formidable inferred (lO'Brien et allboiOt iBaneriee et all 1201 ll) but direct 



X 



; endeavour for some of the most distinguished astronomers and observations of flaring g as disks are technically very difficult, 
mathematicians of the last three centuries to discover the fig- even if DGs are edge-on dSancisi & Allenl 119791) . N evertheless 



ures of equilibrium in the presence of rotation. In particular, mode r ate flaring have b e en ob served in the G alaxy (iKent et al 



the discovery of Jacobi in 1834 that equilibrium figures of uni- |1991l: iKalberla & Kerpl |2009h and in M31 dBrinks & Burton , 



formly rotating fluids need not be axisymmetric took the sci- 1 19841) and it is thus reasonable to expect flaring also in some gas- 

entific community by surprise. The formidable body of knowl- r i c h DGs. Because of the complex geometries (and, often, of the 

edge on (incompressible and uniformly rotating) equilibrium fig- rotation curves) of realistic galaxies, it is extremely complex (if 

ures has been fil l ed up , corrected and consolidated only recently not impossible) to analytically compute figures of equilibrium 

jChandrasekhari 119691: iTassouli [19781). and one must resort to numerical methods. 

Unfortunately, this knowledge has proven to be inadequate An equilibrium model without gas self-gravity suffers from 

for the study of equilibrium configurations in galaxies mainly two major drawbacks. First, such models cannot in principle 

for two reasons: (/) the gas in galaxies is not incompressible, be used to infer equilibrium configurations prone to star for- 

(ii) galaxies are not uniformly rotating. Fortunately enough, mation since the star formation criteria explicitly or implic- 

non-axisymmetric structures of equilibrium are secularly trans- itly rely on self-gravity as one of the key ingredients for star 

formed into axisymmetric figures in realistic (compressible, formation. This lack of self-consistency may lead to situa- 
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tions when the star formation feedback due to supernova ex- 
plosions is studied in models that are insusceptible to star for- 
mation in the first place. Second, neglecting gas self-gravity 
one runs the risk of building a gravitationally overstable con- 
figuration, which would have never been realized if self- 
gravity were taken into account. Such a non-self-gravitating 
configuration would have too much gas compared to the self- 
gravitating counterpart and additional theoretical or empirical 
crit eria are usually invoked to constrain the total ga s mass (see 
e.g. lMac Low & Ferraral Il999l IVasiliev et al.ll2008l) . Moreover, 
the energy rele ase and the corresponding SF rates are often set 
arbitrarily (e.g. iMac Low & Ferraral 11999b . 

This paper is the first of a series of works dealing with the 
dynamical and chemical evolution of gas and stars in DGs em- 
bedded in dark matter (DM) halos. In the context, achieving an 
initial equilibrium configuration is clearly necessary in order to 
study how the onset of an episode of star formation or of an- 
other perturbing phenomenon affects the evolution of the stud- 
ied object. Surprisingly, almost all the papers on this subject ne- 
glect self-gravity of gas and stars and consider a simplified initial 
equilibrium configuration, namely a rotating isothermal gas dis- 
tribution in hydrostatic equilibrium with a fixed potential well (a 



DM h alo or a static distribution of stars; see e.g. Suchkov etaL 
1994 IMac Low & Ferraral Il999t IStrickland & Stevens! 12000 



Recchi et all 12001b iMarcolini et all 12003 ~ IVorobvov et al 
2004 IScannapieco & Briig gen, 2010, among many others). 

In this work, we solve numerically the steady-state momen- 
tum equation of a multi-component galaxy (made of gas, stars 
and a DM halo) taking into account the gravitational accelera- 
tion of all these com ponents. This task has been attempted only 
by very few authors dNaravan & Jog! 12002b lHarfst etafl 12006b 
iBaneriee et al.bl20Tll) . The typical justification forneglecting gas 
self-gravity in constructing initial equilibrium models for DGs 
is that "the gravitational potential o f DGs with M g < IP 9 M Q i s 
dominated by the dark matter halo" (IMac Low & Ferraral [19991 . 
However, it is worth reminding that some authors still doubt 
about the presence of massive DM halos around DGs. For in- 
stance, recent observations of the mass- t o-ligh t ratios in Virgo 
Cluster dwarf ellipticals by Toloba et al.l (1201 ll) and in gas-rich 
DGs bv lSwaters et al.l d201 1 ) and also studies of s tructural prop- 
erties of the Milky Way dwarf spheroidals (see e.g. iKroupa et all 
120101 and references therein) are substantially questioning the 
contribution of DM on small scales. 

Moreover (and more importantly), even if the total mass 
of a dwarf galaxy is dominated by a DM halo, within the 
Holmberg radiu s most of the galaxy i s made of baryons (see e.g. 
iPapaderos et all 1 1996b ISwaters et all 201 lh. although some au - 
thors report different claims (e.g. ICarignan & Beaulieub [19891) . 
In some numerical works (which neglect gas self-gravity) it can 
be clearly noticed that the assumed DM profile leads to a very 
low density of the DM component (much lower than the gas 
density) in the central region of the simulated galaxy (for in- 
stance in iD'Ercole & Brighentib fl999l the central gas density is 
~ 10 times larger than the DM density). Central densities of the 
DM halos in DGs can also be inferred from the observed rota- 
tion curves an d typical values are q uite low; significantly below 
10~ 24 g cnr 3 dde Bloket alll2008l) . We can thus conclude that it 
is very unlikely that gas self-gravity is negligible in the central 
parts of gas-rich DGs. All the above arguments in favour of gas 
self-gravity justify the relevance of the present study. 

The plan of the paper is as follows. The basics of the numer- 
ical model are described in Section [2] The initial and boundary 
conditions, as well as the solution procedure, are summarized 
in Section [3] The main results are presented in Sections |4] and 



[5] A comparison of our results with predictions of the ACDM 
theory is given in Section [6] The implications for the evolution 
of DGs and the model caveats are discussed in Sections|7]and[8] 
respectively. The main conclusions are summarized in Section|9] 

2. Numerical model 

An accurate construction of self-gravitating, rotating equilibria 
involves solving for the steady-state momentum equations of 
gas, stars, and dark matter in their combined gravitational poten- 
tial. This is however a difficult and time consuming numerical 
exercise, since the density distribution of each component de- 
pends on the total gravitational potential, which in turn depends 
on the spatial distribution of each component. We simplify our 
task by making two assumptions. First, we neglect the contribu- 
tion by the stellar component to the total gravitational potential 
throughout most of the paper and return to quantify this effect 
in Section [531 Second, we assume that the DM halo has a fixed 
form and hence a fixed gravitational potential. We note that this 
assumption may break down on timescales much longer than a 
galactic orbital period. We plan to investigate the response of the 
DM halo in a follow-up study. 

The resulting steady-state momentum equation for the gas 
component in the total gravitational potential of gas and dark 
matter takes the following form. 



— VP + (v V)v = g g +g h , 

Pi 



(1) 



where P = p g cr g is the gas pressure, cr g is the one-dimensional 
gas velocity dispersion, p s is the gas volume density, v is the gas 
velocity, and g a and g h are the gravitational accelerations due 
to the gas and DM halo, respectively. The gravitational acceler- 
ation g g is calculated as g g = — VO, where the gas gravitational 
potential <J> is obtained via the solution of the Poisson equation 



V 2 = 4nGp g . 



(2) 



For rotating equilibria, it is most convenient to expand equa- 
tion (LTJ in cylindrical coordinates (r, z) with imposed axial sym- 
metry, i.e., d/d(f> = 0. A steady state solution implies that v z = 
and v t = (but + 0) and the resulting equations are 



7 dlnp g 
Fg dr 

2 dln Pg 

dz 



doi 

~d7 



+ gg,r + gh,r 



do* 
' dz 



+ gg,z + gh,z- 



(3) 
(4) 



In this study, we assume that the gas temperature is spatially uni- 
form (see Section l3~2l i. which implies that the spatial derivatives 
of cr 2 are zeroed. 

Equations (01 and (0]i are discretized using a first-order 
backward-difference scheme on a cylindrical mesh with 600 x 
600 grid points assuming the axial and midplane symmetry 
around the z-axis (r = 0) and the midplane z = 0, respectively. 
The resulting set of 600 2 linear equations is solved using a for- 
ward substitution scheme explained in detail in the Appendix. 

3. Initial conditions 

3. 1 . Dark matter halo setup 

In order to solve equations (0 and (0) for the gas density p g , 
one needs to specify the form of the DM halo. We take two dis- 
tributions that are most often used to fit the rotation curves of 
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DGs. The first choice is a quasi-isothermal sphere, which has 
a flat near-central density distribution and a tail inversely pro- 
portional to the square of the distance from the galactic center 
m - (z 2 + r 2 ) 1 / 2 and is described by the following equation 



Pqis = 



A) 



l + (Wro) 2 ' 



(5) 



The central density po and the characteristic scale length of the 
quasi-isoth ermal halo can be calculated using the following re- 
lation s (e.g. lMac Low & Ferraralll999llSilich & Tenorio-Tagle I 
20010 



x S [AfbM| ffl/2k 



1/2 



po = 6.3 x 10 



10 / ^DM 



-1/3 



H~ l/3 M B kpc~\ 



(6) 



(7) 



where Mum is the mass of the DM halo contained within the 
virial radius 

.oief^VV^kpc. 



(8) 



We note that for a fixed H (set to 0.65 in the current pape r 
for consistency with the work of iMac Low & Ferraral (119991) ) 
the quasi-isothermal halo is uniquely determined by a choice 
of Mqm- Finally, the gravitational acceleration of the quasi- 
isothermal halo can be written as 



g h = -AnGporl/w 2 [vr/r - arctan (ar/r )] e», 



(9) 



where e* = m/wis the unit vector. 

The second choice for the form of the DM halo is the well - 
known NFW density profile suggested by iNavarro et al. I dl997h . 
which features a cuspy profile in the inner regions and a tail in- 
versely proportional to vr" 



Pnfw - 



Pc 



(cr/r c )(l +m/r c ) 2 ' 



(10) 



where p c and r c are free parameters. The mass of the DM halo 
contained within radius m can be expressed as 



M DM (m) 



M, 



DM 



ln(l 



1 + x 



(ID 



m 

where x = wc/nr v i r , f(c) = ln( 1 + c) - c/( 1 + c), c is the concen- 
tration parameter, and the virial radius tu vw is defined by equa- 
tion ([8]l. The concentration parameter is det ermined from the 
statistics of the ACDM halo concentrations bv lNeto etall d2007l) 



c = 4.67 



Mr, 



10 14 M F 



H~ 



(12) 



Finally, the gravitational acceleration due to the NFW halo can 
be calculated as 



GM DM (m) 
gh = o e. 



(13) 



A DM distribution profile, somewhat intermediate between 
the NFW and the quasi -isothermal pro files, has been semi- 
empirically introduced by IBurkertl (1 1995b . It is described by the 
following equation: 



PBurk 



Pc 



(1+Wc)[l + W] 



(14) 



1 We note that IMac Low & Ferraral d 19991) have a misprint in the ir 
equations which has been corrected in lSilich & Tenorio-Tagle I d200lh . 



It is thus a cored profile (as the quasi-isothermal one) which, 
in analogy to the NFW profile, declines at large radii as nr~ 3 . 
Although this profile fits well the rota tion curves of DGs 
(IBurkertl 119951: ISalucci & Burked l2000l) . we have not taken 
it into consideration, because the results adopting this profile 
are intermediate between the results with a (cuspy) NFW and 
a (cored) quasi-isothermal profile. As we show in Sect. 15.21 our 
results depend very little on the DM profile, hence, for the sake 
of conciseness, we have not considered the Burkert profile. 

3.2. The gas temperature 

The thermal properties of gas affect the form of the resulting 
equilibrium configuration. A fully self-consistent approach re- 
quires solving for the thermal balance equation along with the 
steady-state equations © and This however entails a con- 
siderable increase in calculation time and, sometimes, results in 
poor convergence. 

In this study, we take a simpler approach and build equi- 
librium configurations for a pre-defined gas temperature. This 
approach is justified if the characteristic cooling/heating time of 
gas is much shorter than the dynamical time. The pre-defined 
gas temperatures are varied in a wide range, starting from 100 
K, typical for the cold atomic clouds, to a few x 10 4 K, typical 
for the warm diffuse gas. 

3.3. Rotational versus thermal support 

In order to construct rotating equilibria, one needs to specify the 
form of the rotation curve. A usual approach is to set the ro- 
tation velocity of gas v<j to the circular velocity v c j lc , thus as- 
suming th at the support against gravi ty comes mainly from rota- 
tiotFke.g. lMac Low & Ferraral 1 19991) .' Such an assumption taken 
blindly may produce flattened gaseous disks with surface den- 
sities nearly independent of galactic radius or even increasing 
outward, which is unlikely when compared to real systems^. 

A more general and realistic approach is to assume that part 
of the support against gravity comes from pressure gradients and 
to set v<s = ov C irc> where a is the spin parameter that determines 
the relative contribution of rotation to the total support against 
gravity. For or = 1, the gas disk is almost totally supported by 
rotation, whereas for a — the disk is thermally supported. The 
resulting expression for the rotational velocity of gas is 

v = a[r(\ 8h Az = 0)| + \ 8s Az = 0)|)]' /2 , (15) 

where the subscript r denotes the radial component of the gravi- 
tational accelerations, the latter being calculated in the midplane 
z — C0. This choice makes the rotation velocity z -independent, in 
concordance with the Poincare-Wavre theorem (tLebovitzl 1 19671) 
for a barotropic gas in rotation equilibrium. More realistic ro- 
tating equilibria with a negative vertical gradient of requir e 
considering a more general baroclinic gas (iBarnabe et all 120061) . 
which is out of the scope of the present study. Throughout 
most of the paper, we use a = 0.9 dTomisaka & Ikeuchil U988I: 



2 Still there will be some support from gas pressure gradients because 
circular velocity is not an exact solution of the steady-state equation (TJ 
with P + 0. 

3 As mentioned in the Introduction, it is very likely that the gas verti- 
cal scale height naturally increases outward, producing flaring, but the 
vertically integrated gas volume density (i.e., surface density) decreases 
outward. 

4 In practice, we calculate |gh, r (z = 0)| and |g a-r (z = 0)| in the first 
layer of computational cells lying right above the midplane. 
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IStrickland & Stevensl l2000t) and explore the dependence of our 
results on smaller values of a in Section lBTTl 



3.4. Boundary conditions 

The final step is to specify the values of gas volume density p s 
at the boundaries. We use a computational box with physical di- 
mensions of 8.0 kpc along the r- and z-axes, with the spatial 
resolution of 13.3 pc along each coordinate direction. Reflecting 
boundary conditions at the z- and r-axes are a natural choice. 
In addition, one needs to define the bounding pressure, i.e., the 
values of p g and cr g at the outer z- and r-boundaries, if a galaxy 
is submerged in a dense and hot intra-cluster medium. These 
values however are essentially free parameters, since they de- 
pend on the environment. Hence we decide to take a different 
approach and define the value of gas number density «o,o in the 
innermost computational cell near the origin (z = 0, r = 0). 
This value is kept fixed throughout the iterative solution pro- 
cedure (described below) and serves as a "seed" density needed 
to solve equations (0 and (0- Increasing/decreasing the value 
of «o,o would yield more/less massive gase ous disks of different 
spatial configuration. Thus, at variance wifh lMac Low & Ferraral 
(Il999h . our model galaxies do not have a disk cutoff radius. On 
the other hand, because of that, ours are truly equilibrium config- 
urations and the disk does not tend to expand into the intracluster 
medium as in Mac Low & Ferrara (see the test problem in the 
Appendix). 

The choice of a, cr„ and «n,o> along with the mean molecular 
weight fi — 1 .26 (for a metallicity of 1/100 that of the solar) and 
reflecting boundary conditions at the z- and r-axes, completes 
the initial setup and allow us to calculate equilibrium configu- 
rations of gaseous disks for different shapes and masses of the 
DM halo. The bounding effect of the external environment will 
be addressed in a future study. 



SOLUTION SCHEME 



Setting DM halo and gas temperature. 
Initial guess for p init and m 



Calculating 
g g from p init 



Obtaining p from 
eqs (2) and (3) 



Calculating 
g g and from p init 



Convergence 
S < 10- 5 



NO 



Another iteration 
Pinit = P 




Fig. 1. Schematic representation of the iterative solution proce- 
dure to calculate rotating, self-gravitating gaseous equilibria in 
DM halos. 

narrow, high-density plateau which is difficult to resolve numer- 
ically. 



3.5. The solution procedure 

One can notice that, since g g is obtained through the solution of 
the Poisson equation 0, equations and (0 are transcenden- 
tal and hence require an iterative solution procedure described 
schematically in Fig. Q] The calculation begins with a) choos- 
ing the DM halo profile (quasi-isothermal sphere or NFW halo), 

b) calculating the corresponding gravitational acceleration g h , 

c) fixing the thermal properties of gas, i.e., the gas temperature, 
and d) making an initial guess for the gas volume density p g ,i n it 
and rotational velocity v^unt. For the former we usually choose 
a spatially uniform distribution with p g = inn fi «o,o and the lat- 
ter is determined from equation ([TBI . Finally, the gravitational 
acceleration g a of the gas configuration is calculated by solv- 
ing for the Poisson equation using the alternative di r ection im- 
plicit method as describ ed in lBlack & Bodenheimerl d 19751) and 
IStone & Normanl d 19921) . 

The first loop of iterations begins with solving for the steady- 
state equations and 0. The resulting gas density p g is com- 
pared against the initial guess p gj ; n it for every computational cell 
and if the maximum relative error e is larger than 10~ 5 , then 
the iteration cycle repeats by setting p g ,i n it = p g and calculating 
new gravitational potential and rotation curve of the gas disk. 
Usually, convergence is achieved after 10-15 iterations, but may 
sometime diverge signalizing for an inappropriate initial guess 
for «o,o or a, especially when «o,o is large and the corresponding 
equilibrium solution near the rotation axis is characterized by a 



4. Self-gravitating equilibrium gaseous disks 

4.1. Star formation criteria 

We build equilibrium gaseous disks hosted by DM halos of 
various mass and shape and determine the minimum gas mass 
needed to trigger star formation in these systems. Three crite- 
ria are employed to assess the feasibility of star formation in our 
model galaxies. The first criterion is based on theoretical consid- 
erations of gravitational stability in self-gravitating systems. We 
assume that star formation is allowed if the Toomre Q parameter 



ttGE 



(16) 



is smaller than a critical value Q c , where v is the epicycle fre- 
quency and £ is the gas surface density. The classical analysis of 
thin, axisymme tric gaseous disks suggests a value of Q c — 1 .0 
(iToomrel 1 19641) . but Q c in real systems is usually somewhat 
greater and may depend on many factors including the galaxy 
class, the form of the rotation curv e, the disk thickness, the 
strength of magnetic fields, etc (e.g . iPolvachenko et all 1 1997 
| Kim & Ostriken. |200lt iBigiel et all 120081: iLerov etaU 12008 



Dong et all 120081: iRovchowdhurv et all l2009l> . In this study, we 
take a conservative value of Q c = 2.0 and assume that our model 
galaxy is prone to star formation if £>t < Qc in at least some 
parts of the gas disk. 

For the second star formation criterion, we make use of 
empirical studies of star formation in the Local Universe by 
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Table 1. Quasi-isothermal DM halo parameters 





r (kpc) 


po (M Q per 3 ) 


ETvir (kpc) 


10' 


2.27 x 1(T 2 


0.337 


4.6 


10 8 


7.2 x 1(T 2 


0.157 


9.9 


10 9 


0.23 


0.073 


21.3 


10 10 


0.72 


0.034 


45.9 



lKennicuttllll998[l2008h who compares disk-averaged star forma- 
tion (rates per unit area) versus gas surface densities in normal 
and starburst galaxies, including DGs. These studies suggest the 
following scaling law (hereafter, the Kennicutt-Schmidt law) be- 



tween the star formation rate per unit area £sfr(M q yr 1 kpc 2 ) 



and the gas surface density £(M pc 2 ) 



Ssfr = (2.5 + 0.7) x 10" 



1 Mq pc -2 



1.4±0.15 



(17) 



with a threshold density S c of the order of 5.0 M pc 2 , below 
which very rare cases of large-scale star formation are detected. 
Therefore, we use this value as the second star formation crite- 
rion and assume that star formation can be triggered in our model 
galaxies if 2 > E c . 

So far, we have used vertically integrated gas densities to 
assess the model's susceptibility to star formation. However, 
star formation recipes may also rely on the critical gas vol- 
ume density n c , as is often done in num erical hydrodynam - 
ics simulations. Moreover, as discussed in lElmegreenl (ll997h . 
a Schmidt law with index 1.5 would be expected for self- 
gravitating disks, if the SF rate is equal to the ratio of the lo- 
cal gas volume density to the free-fall time, all multiplied by 
some efficiency. The adopted values of n c vary in wide lim- 
its, depending on the numerical resolutio n but most studies use 
value s of the order of 0.1-1.0 cm 3 (e. g. ISpringel & Hernquistt 
l2003t ISchave & Dalla VecchiaL 120081). tho ugh some authors 
adopt much higher values (e.g. lTaskerll20lTl) . 

In this paper, we assume that SF is allowed if there is enough 
gas mass (in the vertical column that fulfils the first two criteria) 
with number density « g greater than a fiducial critical value of 
n c = 1.0 cirT 3 to allow for a SF event of non-negligible magni- 
tude, i.e., if 



M; 



SF 



(»g 



> 1.0 cm 



3 ) > 10 4 



M 



(18) 



We choose to set a limit in mass rather than in size because star 
formation may be localized to just a few tens of parsec, yet con- 
tain enough gas mass for a star formation event of notable mag- 
nitude. Finally, a model galaxy is assumed to be prone to star 
formation only if all three criteria are satisfied altogether. 

4.2. Equilibrium models 

Throughout most of the paper, we use a quasi-isothermal DM 
halo with four different masses M DM = 10 7 Mq, 10 8 M Q , 
10 9 M , and 10 10 M . The corresponding values for ro, po, and 
m are listed in TableQ] We set the gas spin parameter to a = 0.9 
and the gas temperature T g to a spatially uniform value that is 
either independent of the halo mass (T g = 10 4 K) or scales with 

2/3 

the DM mass as T s oc M DM , as suggested by the virial relations. 
We consider the effect of varying rotational support against grav- 
ity (i.e., varying a) in Section 137X1 the effect of a different DM 
halo configuration (i.e., the NFW halo) in Section [5721 and the 
effect of non-negligible stellar disk in Section[ 



To put things in the physical context, the models considered 
here and in Sections 151 45.21 (with a gas distribution in equilib- 
rium with a DM halo, waiting for the onset of star formation) 
can be considered as progenitors of, e.g., Blue Compact Dwarf 
galaxies wh ose stellar populations a re largely dominated by very 
young stars (iPapaderos et al 1 120081) . In Section 15731 we describe 
objects with a pre-existing disk of stars which have smoothly 
accreted gas and have achieved a new equilibrium configuration 
(galaxies surr ounded by extended g as reservoirs are quite com- 
mon (see e.g. Ivan Zee et all 1 1 9981 for the case of I Zw 18)). 
Finally, in Section 15741 we assume a redshift of galaxy formation 
significantly larger than zero. Therefore, our equilibrium con- 
figurations should be treated as proxies to DGs that have built 
up their gas mass reservoir by a quasi-steady accretion or have 
temporally achieved a quasi-steady state after an episode of fast 
accretion. 

Figure [2] presents gas surface densities £ (left column), 
Toomre Q parameters (middle column), and gas rotation veloci- 
ties i'0 (right column) for various steady-state gaseous disks with 
Mdm ranging from 10 M (top row) to 10 7 M (bottom row). 
The spin parameter and the spatially constant gas temperature 
are the same for all models and are equal to 0.9 and 10 4 K, re- 
spectively. In the calculation of the Q parameter and we use 
mass-weighting according to the gas mass contained in every 
computational cell. The gas surface density is obtained by inte- 
grating p g along the z-axis. 

For each value of Mdm, we consider five models with dif- 
ferent seed values of the gas number density «o,o, namely, 0.01, 
0.1, 1.0, 5.0, and 25 cm 3 . These models are distinguished in 
Figure [2] by lines of different style, with the dotted lines corre- 
sponding to «o,o = 0.01 cirT 3 and dashed lines to «o,o = 25 cm -3 
(and the other models in between in the order of increasing «o,o)- 
The largest/smallest values of «o,o produce models with high- 
est/lowest gas surface densities near the galactic center. 

The radial profiles of £ in Figure [2] indicate that models 
with lower values of Mdm produce less centrally concentrated 
gaseous distributions. Indeed, models with Mdm = 10 10 M Q 
have a density tail proportional to r 2 , whereas models with 
Mdm < 10 8 M are characterized by S oc r _u . This tendency 
can be explained by the fact that the mass of the gas disk starts to 
systematically exceed that of the DM halo for M D m < 10 9 M , 
the effect discussed in more detail in Section [5] As a result, the 
shape of the gas disk in models with Mdm < 10 9 M is mostly 
determined by self-gravity of the gas, with the resulting distri- 
bution approaching that of a self-gravitating isothermal ellipsoid 
with the density tail p g oc m~ 2 or X oc r . We also note that mod- 
els with lowest values of X tend to have surface density profiles 
independent of radius. 

For a given DM halo mass and equal gas temperature, mod- 
els with lower values of £ have higher values of the Q parameter, 
as expected. It is worth noting that the lowest Qj is often found 
a few hundred or even thousand parsecs away from the galactic 
center. This behaviour can be understood by analysing the radial 
dependence of the epicycle frequency v = (4Q 2 + rdQ. 2 /dr) 1 ^ 2 . 
This quantity is independent of radius r in the inner parts r «: ro, 
where the DM halo and gas densities are nearly constant and 
Q. ss const. On the other hand, at r » ro the epicycle frequency 
declines with radius because the DM halo and gas densities also 
(as a rule) decline with radius. This implies that Qj is nearly in- 
dependent of r in the inner parts but may increase or decrease in 
the outer parts depending on the radial profile of the gas surface 
density E. For models with £ nearly independent of radius, Qt 
generally decreases at large radii (because v also decreases but 
other quantities stay nearly constant), whereas for models show- 
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Fig. 2. Gas surface densities (left column), Toomre Q parameters (middle column), and rotation velocities (right column) of self- 
gravitating (steady-state) gaseous disk hosted by quasi-isothermal DM halos with four masses Mqm as indicated in the right column. 
The spatially uniform gas temperature and the spin parameter are set to T g = 10 4 K and a = 0.9, respectively. The horizontal dotted 
lines mark the adopted critical surface density for star formation S c — 5 M Q pc 2 (left column) and the critical Toomre parameter Q a 
for gravitational stability (middle column). For every value of Mum, five models with different values of the seed density «o,o=25.0, 
5.0, 1.0, 0.1, and 0.01 crrT 3 from the uppermost to the lowermost line are considered. Greater «o,o produce models with higher gas 
surface densities in the inner regions. The models that are susceptible to star formation according to the adopted star formation 
criteria are marked with red thick lines and the gas mass M^ F with number density greater than a critical value of 1.0 cm~ 3 is 
indicated in the left column for models plotted with solid red lines. 



ing a decline in S at large radii, the corresponding values of Qj 
attain a minimum at some several hundred or thousand parsecs 
and increases on both sides. 



It is also worth noting that the rotation curves of our model 
galaxies either steadily rise or flatten out only at large radii. This 
behaviour is in qualitative agreement w i th the observed rota- 
tion curves of DGs (see e.g. lde Blok et aTl l2008tlOh et"aTll2008t 
ISwaters et all I2009L 1201 lh . 



The horizontal dotted lines in the left and middle columns 
of Fig. |2]mark the critical gas surface density for star formation 
E c = 5 M e pc 2 and the critical Toomre parameter Q c = 2.0. We 
use these values to determine models that can allow for star for- 
mation, i.e., models for which both criteria X > E c and Qj < Q c 
are met at least in some parts of the galactic disk, and also the gas 
mass F with number density greater than 1 .0 cirT 3 exceeds a 

minimum value of 10 4 M as stipulated by the third SF condi- 
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Fig. 3. Total gas masses M g calculated for five models with 
different seed density values «o,o- The filled squares mark the 
models that are susceptible to star formation according to the 
adopted criteria. Lines of different style connect models with 
equal masses of the DM halo. 



tion. Such "star-formation-allowed" models are highlighted with 
red color in Fig. [2] 

For every model in Figure |2] we calculate the total gas mass 
M„ contained within our computational domain, the latter having 
a cylindrical shape with radius of 8 kpc and height of 8 kpc on 
both sides from the midplane. We use M g as a proxy for the 
total gas masfl and estimate the possible missing gas mass using 
the surface density profiles in Fig. [2] For models with Mdm = 
10 10 M Q , £(r) oc r~ 2 and M g (r) oc ln(r) at large radii, implying 
a small correction of order unity for a computational box with 
size three times greater than ours (24 x 24 kpc). For models with 
M DM < 10 s M Q , Z(r) oc r~ u and M g (r) oc r 09 , implying a 
factor of 2.7 increase in the total gas mass. This means that our 
estimates are accurate to within a factor of unity for models with 
Mdm ~ 10 9 M , while for models with smaller DM halos we 
may underestimate the total gas mass by up to a factor of 3. 

The total gas masses M g for every model are plotted in Fig. [3] 
The bottom axis shows the seed number density «n,o f° r eacn 
model. Lines of different style connect models with the same 
mass of the DM halo (e.g., models with Mdm = 10 9 M Q are 
connected with the dash-dotted line). We distinguish the star- 
formation-allowed models by the filled squares. 

It is evident that DM halos can host steady-state gaseous 
disks with various masses, but not all gas configurations are 
prone to star formation. There exists a minimum gas mass 
M™ 11 that a DM halo needs to accumulate in order to fulfil the 
star formation criteria. For instance, for a DM halo with mass 
10 10 M (dotted line) the corresponding minimum gas mass is 
4.3 x 10 8 M , while for a DM halo with mass 10 8 M Q (dashed 
line) the corresponding minimum gas mass is 1.3 x 10 9 M . As 
one can see, M™ n depends on the mass of the DM halo and may 
increase as Mdm decreases. The latter effect is not unexpected 
- a DM halo with smaller mass has a weaker gravitational po- 
tential and, as a consequence, a more massive gaseous disk is 
required to attain the critical density for star formation. The gas 

3 We cannot extend our computational boundaries to a distance much 
greater than 8 kpc because of the need for high spatial resolution to 
resolve the inner density plateau. 



self-gravity here is a key factor, without which such an effect 
will be absent. 

5. Minimum gas mass for star formation vs. dark 
matter halo mass 

In this section, we study in more detail the dependence of the 
minimum gas mass for star formation M™ n on the mass of the 
DM halo, as well as on other properties of galactic systems. 
We want to emphasize here that M™ n is the total gas mass of 
a galaxy in which star formation is allowed according to the 
adopted SF criteria and not the gas mass M 8F that fulfils the star 
formation criterion (IT8V The latter quantity is always smaller 
than M™ m as only part of the gas disk is characterized by n g > « c . 
We choose to concentrate on M™ n because we compare this 
quantity to the baryonic mass derived from the ACDM theory. 

Throughout the paper we consider quasi-isothermal DM ha- 
los (if not specified otherwise) described by equations (|5}-©. 
The top panel in Figure [4] presents the M™ n -MoM relation for 
the four values of the DM halo mass (10 7 M , 10 8 M Q , 10 9 M , 
and 10 10 M ). In particular, the thick blue solid line shows 
the data for the spin parameter a = 0.9 and gas temperature 
T g = 10 4 K (independent of the DM halo mass), whereas the 

2/3 

thick blue dashed line does that for T g oc M DM and the same 
value of 00. The latter relation is normalized to T g = 10 4 K for 
Mdm = 10 10 M , which yields the following scaling law 

r g =2.15xl(r 3 M D M- (19) 

The corresponding values of r a are indicated in Fig. |4]for 
every pair of data (M™",Mdm)- The adopted set of parameters 
(quasi-isothermal DM halo, a = 0.9, and gas temperature either 
dependent on or independent of Mdm) is denoted hereafter as 
the reference model. The bottom panel in Figure |4] also shows 
the ratio £ = M™"/M D m versus M D m- 

The dotted line in Figure |4] presents the baryonic mass Mb 
for a given DM hal o mass as exp e cted f rom the ACDM theory 
and WMAP4 data dSpergel et all 120031) with Q b /Q m = 0.17. 
If we treat Mb as an upper limit of the available gas mass, it 
becomes evident that galaxies with Mdm < 10 9 and T g = 10 4 K 
require much more gas than available to achieve a state in which 
star formation is allowed. This statement applies to quasi-steady 
systems and may break down for galaxies that accumulate their 
gas mass reservoir via a series of violent mergers or if an external 
perturbation drives the system out of equilibrium and triggers 
star formation in some parts of the galaxy, as discussed later in 
Section O 

An alternative solution is that galaxies can cool to suffi- 
ciently low temperatures to warrant a more compact and dense 
gas configuration. Figure [5] shows gas surface densities £ (left 
column), Toomre Q parameters (middle column), and gas rota- 
tion velocities (right column) for the same parameters as in 

Figure [2] but for T g oc as described by eq. $1% , with the 

corresponding gas temperatures indicated for every Mdm in the 
right column. One can see that the gas surface density profiles 
are considerably steeper for cooler gas disks and are character- 
ized by approximately the same power law r~ 2 in the outer re- 
gions. Furthermore, the transition radius from a near-constant 
surface density to the sloped one decreases with masfl in con- 

6 We remind that in all models the gas temperature is spatially uni- 
form but its value may or may not depend on Mdm- 

7 Diminishing pressure support against gravity in models with 
smaller T e is compensated by an increase in the gas density slope. 
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Fig. 4. Top. Minimum gas mass for star formation M™ n as a 
function of the DM halo mass Mdm- The thick blue lines show 
the data for the reference model (quasi-isothermal DM halo, spa- 
tially uniform T g , and a = 0.9), while the thin red lines show the 
effect of a smaller spin parameter a = 0.5. The solid lines (of 
corresponding color) refer to models with T g = 10 4 K, while 
the dashed lines show models with gas temperature according 
to eq. ( fT9l . Numbers denote the corresponding gas temperatures 
for every model. The dotted line gives the available baryonic 
mass Mb as suggested by the ACDM theory with Qb/Q m = 0.17. 
Bottom. Ratio M™ n /M DM versus M DM for the same models as 
in the top panel. 



trast to models with the DM-mass-independent gas temperature 
Tg (Fig. 2). 

The star formation criteria (ITol and (TlTt are satisfied only in 
the inner parts of our model galaxies, with the size of the star 
formation region shrinking to a few tens of parsecs for models 
with M DM < 10 8 M . Moreover, the value of M g F notably de- 
creases for lower mass DM halos. These very compact starburst 
regions are not unusual in low mass star forming galaxies. For 
instance, the galaxy S BS 0335-052 has a st ar forming radius of 
only ~ 20 pc (see e.g. lTakeuchi et all l2003h . 

As the blue dashed line in the top panel of Fig. |4]demon- 

2/3 

strates, models with T g oc M DM have much smaller M™ 11 than 

2/3 



models with T g - 10 K. Moreover, models with T g oc M DM 
are characterized by M™ n < Mb, meaning that such galaxies 
may have enough gas reservoir to achieve a state in which SF is 
allowed. However, the required gas temperatures are quite low, 
especially for Mdm ^ 10 8 M Q . Moreover, the obtained rotation 
curves (third column in Figure [5} are much flatter for models 
with Mdm ^ 10 8 M Q than for their more massive counterparts. 



5.1. The effect of varying rotation support 

Decreasing the amount of rotation leads to more compact and 
dense gas configurations as the pressure forces start to play an 
ever increasing role in the support against gravity. Therefore, 
one can expect that lowering the spin parameter would produce 
steady-state gaseous configurations with steeper gas surface den- 
sity profiles (provided that the gas temperature is constant). 

This effect is illustrated in Figure [6] showing the gas surface 
density distributions for models with a = 0.5 (solid lines) and 
a = 0.9 (dashed lines). For the sake of convenience, we com- 
pare only the star-formation-allowed models^ with the smallest 
value of M™ n . In the case of a — 0.9, such models are distin- 
guished by red solid lines in Figures [2] and As one can see, 
models with smaller a are characterized by more compact gas 
disks with steeper density profiles at large radii. This effect takes 
place because the a = 0.5 models compensate for a smaller ro- 
tation support with steeper gas density (and pressure) gradients. 
As a result, these models also have smaller M g F than the corre- 
sponding a = 0.9 modelfl 

Figure |4] illustrates the effect of a smaller rotation support 
against gravity. The thin red lines show the data for a = 0.5, 
with other parameters being identical to those of the reference 
model. It is evident that the minimum gas mass for star formation 
M™ n is substantially lower in galaxies with smaller values of a. 

The effect is particularly strong for galaxies with T g oc Mj^. 
Lowering further the spin parameter to a — 0.25 yields only a 
factor of two decrease in M™ n at most. We point out that the 

minimum gas mass for star formation M™ n remains still much 
greater than the available baryonic mass Mb for galaxies with 
M DM < 10 8 Mn and T e = 10 4 K. 



5.2. Variations in the DM halo form 



As already mentioned, the mass and shape of DM halos in DGs 
are still very uncertain and the problem is highly debated in the 
literature. If the DM density at small galactocentric radii m is 
approximated by a power law (pdm K et^X a value of /? close to 
or slightly smaller than one (i.e., cusps) is obtained in numerical 
simul ations based on the ACDM theory (iNavarro et al l fl997l 
I2010h . On the other hand, observations suggest values of /3 close 
to zero, implying t he presence of cores with near-constant DM 
density at small m (Ide Blok et all l2Q01b Ispekkens et all 120051) . 
This unsolved mismatch between observations and models is 
commonly referred to as the "cusp-core problem". 

Our reference quasi-isothermal DM halos are cored DM pro- 
files. In this section we estimate the effect that a cuspy DM 



8 Strictly speaking, the M DM = 10 7 M Q and T g = 10 2 K case has 
no star-formation-allowed models since the third criterion, i.e., M| F > 
10 4 M is not fulfilled even for the highest n 00 model. We relax this 
requirement in this particular case since the other two SF criteria are 
nevertheless satisfied. 

9 We note that some of the a = 0.5 models have higher gas densities 
in the innermost regions, which is the result of a rather coarse grid of 
models considered in the paper. Nevertheless, the size of this region is 
rather small (~ 0.2 kpc) and most of the gas mass is still concentrated 
at large radii. 
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halo advocated bv lNavarro et al.l (l!997l) and described by equa- 
tions (fT0b-(fT2l> may have on the values of M™ n . Figure |7]com- 

pares the minimum gas masses M™ n obtained in the reference 
model for the quasi-isothermal DM halo (blue lines) with those 
calculated for the NFW halo (red lines). As one can see, the dif- 
ference is negligible for models with T g = 10 4 K and is minimal 
for models with T g oc M™. Therefore, the form of the DM halo 
is not expected to play a significant role in determining the min- 
imum gas mass for star formation as far as the total mass of the 
DM halo is the same. 

5.3. The effect of pre-existing stellar disk 

So far we have considered model galaxies that consist of a 
gaseous disk and DM halo. However, real star-forming galaxies 
in the Local Universe almost always have a pre-existing stellar 
disk, which may affect the form of the gaseous disk via the stellar 



gravitational potential. To explore the extent of this effect, we as- 
sume that our galaxies in the reference model have a burst of star 
formation that turns 10% of the gas content into stars. Star for- 
mation takes place in parts of the gas disk that obey the three cri- 
teria for star formation Z > Z c = 5 M pc~ 2 and Qj < Q CI = 2.0. 
We then construct a new equilibrium gas disk in the combined 
gravitational potential of gas, stars, and DM halo. The spin pa- 
rameter of the stellar disk is assumed to be equal to that of the 
gaseous disk. This assumption is justified if the gas sound speed 
is comparable to the stellar velocity dispersion, which is often 
true for DGs. 

The black lines in Figure [7] show the minimum gas mass 
M™ n for new models that take into account a recent burst of 
star formation. It is evident that the stellar disk has little ef- 
fect on the value of M™ n irrespective of the DM halo mass (the 
black line is almost indistinguishable from the blue line show- 
ing the reference model). This is explained by the fact that the 
mass of the stellar disk M s is only a small addition to the total 
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Fig. 6. Comparison of the gas surface density profiles for models 
with the spin parameter a = 0.5 (solid lines) and a = 0.9 (dashed 
lines). The DM halo mass Mqm and the gas temperature T g are 
indicated in every panel. The gas mass F with number density 
greater than 1.0 cirT 3 is given only for the a = 0.5 models. 



mass budget. For instance, in models with Mqm = 10 10 M Q , 
M s = 0.07 Mf R and M s = 0.0028 M DM - In models with 
M DM = 10 7 M Q ° M s = 0.047 Mf° and M s = 4.1 M DM (in both 

cases, T g = 10 4 K). Unless many repetitive bursts of star for- 
mation take place with the integrated star formation efficiency 
considerably exceeding 0. 1 (leading to a significant increase in 
the star to gas mass ratio), we do not expect the gravitational 
potential of the stellar disk to affect our results. However, we 
should note that the stellar feedback may drive the gas disk out 
of equilibrium, thus affecting our estimates of M™ n . 



5.4. Galaxies at higher redshifts 

In this section, we study the dependence of M™ n on the redshift 
of galaxy formation z g f. DM halos of the same mass at larger 
redshifts are more compact and one may expect that this could 
affect the shape of the gas disk and hence the value of M™ n . We 
modify equations ©-((S]) to include the dependence on z g t (e.g. 
iFujitaet all 120031) 
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Fig. 7. Top. Minimum gas mass for star formation M™ n ver- 
sus DM halo mass Mqm for the reference model (blue lines), 
models with a cuspy DM halo (red lines), models with pre- 
existing stellar disk (black lines), and models at redshift z g f = 2.0 
(green lines). The solid lines show models with the gas tem- 
perature T g = 10 4 K, whereas dashed lines do that for models 

with T g oc (see equation [19J- Numbers indicate the corre- 
sponding gas temperatures for every model. The dotted line plots 
the cosmological relation between the baryonic and DM mass 
as suggested by the ACDM theory. Bottom. The corresponding 
M™ n /M DM versus M DM relations. 
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where On = "(Zaf = 0) and Cl(z g f) is defined as 
",«(! +z gf ) 3 



"(Zgf) = 



",«(! +z gf ) 3 +Q A ' 



(22) 



(23) 



with Cl,„ and Qa set to 0.24 and 0.76, respectively. 

The green lines in the upper panel of Figure |7]show the min- 
imum gas mass for star formation M™ m as a function of Mum in 
the reference model but with z g f = 2.0. In compliance with the 
downsizing, low-mass galaxies have a median redshift of star 
formation smaller than large objects and a value of z g t = 2 well 
represents an average galaxy formation reds hift for DGs (see e.g. 
iThomas et all l2005t Icattaneo etaTl l201ll) . We also note that 
Zgf = 2.0 yields roughly a factor of two more compact and a 
factor of seven denser DM halos. 

As one can see, the effect of a higher redshift is quite pro- 
nounced for models with T g oc M~ ' (shown by the green dashed 
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lines), producing almost a factor of 10 lower values of M™ n than 
in the reference model (blue dashed lines). However, for models 
with T g — 10 4 K, the resulting values of M' g m quickly approach 
those of the z g f = models (blue lines) for low values of Mqm- 
The lowering of M™ n can be understood if we compare the 
radial density profiles of the DM halos of equal mass at differ- 
ent redshifts. The right panel in Figure [8] presents the volume 
density p q ; s of the quasi-isothermal DM halo as a function of ra- 
dial distance m for the reference model (z g t = 0, dashed lines) 
and the z g f = 2.0 model (solid lines). The corresponding masses 
of the DM halos are indicated in each panel. One can see that 
the DM halos at z g f = 2.0 are characterized by higher densi- 
ties than their z g f = counterparts. This causes gaseous disks 
at higher redshifts to assume steeper density profiles to compen- 
sate an increasing gravity force of the DM halo. This effect is 
illustrated in the left panel of Figure [8] which compares the gas 
surface density profiles in the star-formation-allowed models at 
Zgf = (dashed lines) and the Zgf = 2.0 model (solid lines). All 
other model parameters are those of the reference model and the 
gas temperature is indicated in each panel. Indeed, the z g t = 2.0 
models are characterized by significantly more compact gaseous 
disks than their z s f = counterparts, which results in systemati- 
cally lower values of M™ n . 

6. Reconciling the gas to DM ratio with the ACDM 
predictions 

Figures|4]and|7]demonstrate that the minimum gas mass for star 
formation M™ n may greatly exceed the DM halo mass (and, 
of course, the gas to DM mass ratio can be largely above the 
suggested cosmological value of 0.17) for models with Mum ^ 
10 9 M Q and T g = 10 4 K. This apparent contradiction between 
the baryon-to-DM ratios deduced by our models and the cosmo- 
logically inferred one should not be a major source of concern 
at the present stage. In fact, we know already that some objects 
in the Universe do have baryon-to-DM ratios larger than 0.17 
(e.g., globular clusters, high velocity clouds, and probably also 
dEs and dlrrs, see the Introduction). At the same time, Figures|4] 
and |7] show that models in which T g oc (as suggested by 

the virial relations) neatly reproduce a constant gas to DM mass 
ratio, although below the cosmological value of 0. 17. 

In this section, we explore whether the cosmological value 
can be reproduced by our models with either a different choice 
of the zero point of the T g -Muwi relation or by introducing some 
non-thermal support against gravity in equation ([1}. To explore 
the first possibility, we choose the following scaling law between 
the gas temperature and the DM halo mass 

r g = 4.8 x 10- 3 M^, (24) 

which yields roughly a factor of two higher gas temperatures 
than equation ([19), i.e., T g = 2.23 x 10 4 K for M DM = 10 10 M Q , 
7g = 0.48 x 10 4 K for M DM = 10 9 M , T g = 1.03 x 10 3 K for 
M DM = 10 8 M , and T g = 223 K for M DM = 10 7 M Q . 

The top panel in Figure [9] shows the minimum gas mass for 
star formation M g m as a function of the DM halo mass Mom- 
The dashed line presents the data for the new scaling law de- 
scribed by equation (l24l . while the solid line does that for the 
old scaling given by equation (|19t . As one can see, the new scal- 
ing law yields somewhat higher gas masses than a cosmological 
value and the slope of the model M™-Mdm relation is some- 
what shallower than the cosmological one (dotted line). On the 
other hand, the old scaling law yields M™ n that are somewhat 
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Fig. 8. Left column. Comparison of the gas surface density in 
the reference model (dashed lines) with that of the z g f = 2.0 
model (solid lines). The dotted line marks the critical surface 
density for star formation S c = 5 M e pc~ 2 . Gas temperature 
is indicated in every panel and the DM halo mass is shown in 
the corresponding panels on the right. The gas mass M g F with 
number density greater than 1.0 cirT 3 is given only for the z g f = 
2.0 models. Right column. Volume density profiles of the quasi- 
isothermal DM halos of various mass (as indicated in each panel) 
at Zgf = (dashed lines) and z g t = 2.0 (solid lines). 



smaller than the cosmological values. This means that varying 
the zero point of the Tom k ^dm ration, one can in principle 
achieve a good agreement with the ACDM theory. 

To explore the second possibility, we assume that our model 
gas disk has a non-thermal support against gravity in the form 
of turbulence and magnetic pressure. For the sake of simplicity, 
we assume equipartition between the gas kinetic p ressu r e P an d 
the two sources of non-thermal support (but see ICoxl (120051) ). 
which yields the effective gas pressure P e ff = 3P = 3p g <x g 
(which would correspond to an effective velocity dispersion 
c e ff = y[PesTpl - V3o" g ). The dash-dotted line in the bottom 
panel of Figure [9] presents the M™ n -MoM relation for the case 
with the non-thermal support. It is evident that the correspond- 
ing gas masses increase considerably, in particular for models 
with Mdm ^ 10 9 Mq. The matter is that introducing the non- 
thermal support we effectively increase the gas pressure and the 
corresponding gas surface density profiles become shallower as 
compared to those without non-thermal support. This results in 
an overall increase in the gas mass of a steady-state gaseous con- 
figuration. Figure [9] also suggests that varying the magnitude of 
the non-thermal support one can achieve a good agreement with 
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Fig. 9. Top. Effect of changing the zero-point of the T g oc M DM 
scaling law. The solid lines present the M™ n -MoM relation for 
the reference model with scaling described by equation ([19}, 
whereas the dashed line does that for the scaling given by equa- 
tion (l24l) . The corresponding gas temperatures are indicated for 
each model. Bottom. The effect of non-thermal gas support 
against gravity. Solid lines show the M™-Mdm relation for the 
reference model (zero non-thermal support), whereas the dashed 
line does that for the models with the effective gas pressure three 
times greater than the gas kinetic pressure. The dotted line plots 
the cosmological relation between the baryonic and DM mass as 
suggested by the ACDM theory. 
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Fig. 10. Comparison of the gas volume density profiles along 
the r-axis (left column) and z-axis (right column). In particular, 
models in the top row are characterized by different gas tem- 
peratures T g (as indicated), the second top row by different spin 
parameters a, the third top row by different effective gas disper- 
sions <x eff , and the bottom row by different redshift for galaxy 
formation z g f . All other parameters for every model are the same 
(see text for more details). 



the cosmological mass of baryons given by the dotted line, par- 
ticularly for galaxies with Mqm ^ 10 10 M Q . 

7. Implications for the evolution of dwarf galaxies 

In this work, we have presented numerical solutions for equi- 
librium configurations of model galaxies made up of gas, stars 
and a DM halo in the combined gravitational potential of each 
of these components. The properties of these equilibrium con- 
figurations and, in particular, the minimum gas mass needed to 
achieve a state with allowed star formation M™ n , have been con- 
sidered in detail. 

In future more detailed works, we will be using our derived 
equilibrium configurations as initial setups of galaxy models for 
which we will numerically study the detailed chemical and dy- 
namical evolution, and also the effect of SF feedback. The in- 
terest in simulating the evolution of DGs is steadily growing 
in the last years. The reason is that ACDM theories of struc- 
ture formation predict that dwarf galaxy-sized objects are the 
first virialized structures in the Universe. Moreover, the study 
of star formation and feedback in DGs is in many respects 
much simpler than in large spiral galaxies. Although studies 



of DGs in a cosmological context are more numerous and de- 
tailed than in the past (e.g. iKazantzidis et all 1201 lt | Sales et al.L 
1201 ll iPilkington et all 1201 it iGovernato et all 1201 0l) . still thev 
do not have enough spatial resolution to analyze in detail the in- 
ternal evolutionary processes. A lot in resolution can be gained 
by zooming in and re-simulating small chunks of a l arge cos- 
mological box dMartig et all 120091: ISawala et all 1201 lb . but still 
the best way to accurately simulate a dwa rf galaxy is by numeri- 
cally studying it as a sing l e isolated entity (|Schroven et a"Q 1201 1[ 
IScannapieco & Briiggenl 120101: iRevaz et aL L 120091) although in 
reality they are subject to various environmental effects. 

Although a quantitative comparison between our predictions 
and observations in DGs requires taking into account the feed- 
back from star formation, yet our modeling can give us some in- 
sight as to the expected evolution of DGs. First, we note that for 
some models (especially the low-mass ones with T=10 4 K, see 
Figure|2] bottom row of panels) star formation is not expected to 
occur in the center of the galaxy, but only in a shell with inner 
radius between ^ 100 pc and ^1.0 kpc. 

Second, the evolution of SN-driven shells is known to de- 
pend on the gas density distribution, which in turn is sensitive 
to the initial conditions in a dwarf galaxy. For instance, the 
Rayleigh-Taylor instability in the shell grows faster for steeper 
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gas density profiles. Figure [10]compares the gas volume density 
distributions p g in star-formation-allowed models with different 
T g , a, cr e ff, and z g f. In particular, the left/right columns show 
the radial profiles of p g taken along the r-/z-axes. For the sake 
of conciseness, we consider only models with Mdm = 10 9 M Q , 
models with other DM halo masses show a similar behaviour. 
Models in the top row are characterized by a = 0.9, cr ei j = 0, 
and z g f = (but different T g as indicated in the Figure), models 
in the second top row by T g = 10 4 K, <x e ff = 0, and z a f = (but 
different a), models in the third row by T g — 10 4 K, a — 0.9, 
and Zgf = (but different cr ef f), and models in the bottom row by 
T g = 2150 K, or = 0.9, <x eff = (but different z gf ). 

It is evident that taking a smaller gas temperature T g or spin 
parameter results in equilibrium gas disks with a steeper tail of 
the gas volume density distribution. A similar effect is observed 
for models with a higher redshift for galaxy formation z g f . On 
the other hand, models with and without non-thermal support 
have similar radial profiles of p g . For every model considered, 
the vertical profiles of p g are steeper than those taken along the 
horizontal axis, suggesting a blow-out effect along the rotational 
axis. The variety of possible model realizations implies that the 
evolution of DGs after the onset of star formation may follow 
different pathways depending on the initial conditions in the gas 
disk, even for the same DM halo mass. 

8. Model caveats 

8.1. Steady-state gaseous disks 

A steady-state model is a first-order approximation to DGs 
galaxies. Various effect such as stellar feedbacks, non- 
axysimmetric density waves, and, in particular, external pertur- 
bations may drive DGs out of equilibrium. These phenomena 
can trigger star formation in otherwise quiescent gas disks and 
affect our derived values of M™ n . In order to estimate the possi- 
ble magnitude of such effects, we focus on perturbations with the 
conservation of the total gas mass[3 and refer to star-formation- 
inactive models plotted in Figure |2]by black lines. It is evident 
that positive perturbations in S by a factor of 5-100 are needed 
to drive these models to the star formation threshold. 

If perturbations of such amplitude are possible, then the crit- 
ical gas mass required for star formation may be significantly 
lower. Indeed, as Figure[3]demonstrates, the total gas mass M g of 
an equilibrium configuration declines with decreasing «o,o- The 
filled squares mark the star-formation-allowed models, while the 
open circles correspond to the star-formation-inactive ones. If 
the «o,o = 0.0 lcrrT 3 models can be pushed to the star forma- 
tion threshold, then the minimum mass for star formation M™ n 

may be as low as 10 7 M for the M DM = 10 10 M Q model. 
This corresponds to almost a factor of 40 decrease in the value 
of M™ n as compared to the star-formation-allowed model with 

«o,o = 5 cm 4 . We note, however, this effect becomes con- 
siderably less pronounced for models with smaller DM halo 
masses. For instance, the corresponding decrease in M™ n for the 

Mdm = 10 7 M Q model is only a factor of 2.5. We conclude that 
high-amplitude density perturbations of the equilibrium state can 
significantly affect our estimates of M™ n only for models with 
M DM £ 10 9 M Q but are of rather low significance for models 
with M DM £ 10 s M . 



Perturbations without conservation of the total mass, such as merg- 
ers or ram pressure stripping, can obviously affect our conclusions by 
changing the total available gas mass budget. 



8.2. Multi-phase interstellar medium 

In this study we have neglected the fact that the interstellar 
medium consists of various phases with usually different temper- 
atures and considered a single-phase medium with some typical 
temperature T s . Although cores of molecular clouds (where star 
formation occurs) have temperatures much lower than 10 4 K, the 
latter value must be seen as a mass-weighted mean temperature 
within each computational cell (which has a size much larger 
than the cores of molecular clouds). Indeed, the hot (~ 10 6 K) 
and cold (~ a few x 10 K) gas phases usually amount to about 
1% and 10% of the total gas mass reservoir, respectively (for 
the SF efficiency of ~ 10%), meaning that the mean tempera- 
ture is mostly determined by the warm gas phase. In this sense, 
T g = 10 4 K represents a sort of an upper limit because it is im- 
possible that a computational cell hosting (star-forming) regions 
with temperatures of few tens of K can have an average temper- 
ature significantly above 10 4 K. We note that lower than 10 4 K 
mean temperatures are, of course, possible, provided efficient 
cooling and low SF feedback. 

One may argue that even though some models for T g = 
10 4 K may have difficulty to achieve critical densities for star 
formation, the differentiation into a multi-phase medium may 
eventually push a local region to star formation and that may 
trigger a chain reaction through the bulk of a galaxy. To account 
for this possibili ty, we have adop ted an empirical star forma- 
tion threshold by iKennicuttl (12008 ) . While the Toomre Q param- 
eter criterion is based on the gravitational properties of a single- 
phase medium, the Kennicutt's criterion is based on observations 
of real multi-phase galaxies and hence takes implicitly into ac- 
count the possibility of p hase differentiat ion discussed above. In 
a comprehensive review, iHenslerl d2008l) has discussed the ad- 
vantages of a multi-phase treatment of the ISM in star-forming 
galaxy disks and emphasized the limitations of single gas-phase 
description. 

8.3. Star formation criteria 

In this paper, we have adopted three SF criteria, which are based 
on theoretical arguments, i.e., the Toomre gravitational stabil- 
ity criterion ([Tol , and empirical evidence, i.e., the Kennicutt- 
Schmidt law (fTTt . These SF criteria are not without limita- 
tions and observations suggest that SF can occur even in the 
Toomre-stable regions with Qj ^ Q c , e.g., near the galactic cen- 
ter where the gravitational stability may be determined by the 
rate of sheer rath er than by the magnitude of epicycle motions 
(IVorobvovi 12001 . Moreover, a few galaxies in the Kennicutt's 
(2008) sample harbour SF below the adopted gas density thresh- 
old of E c =5 M pc~ 2 . At the lower end of the KS correlation 
the scatter of measurements also widens, because the SF rate 
fluctuates more stochastically and e.g. starbursti ng DGs are sys - 
tematically located above the relation (see e.g. IHenslerl |2012|) . 
In addition, it is well documented that the KS relation is tighter 
when the SF rate is correlated with the m olecular hydrogen H2 
dKennicutt et aUl2007tlBTgiel et alj . l2008h . 

Nevertheless, we can use the same arguments as in 
Section [8Tl to show that a factor of ten variation in the adopted 
value of E c can significantly affect our estimates of M™ n only in 

models with the DM halo mass M D m ^ 10 9 M Q and the effect 
of uncertainty in E c is diminishing for DM halos with smaller 
mass^ The value of n c is more uncertain and depends largely 
on numerical resolution. Our adopted value of 1.0 cm -3 com- 

11 Note that g c allows for a smaller variation of order unity. 
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plies with most numerical studies on galactic star formation and 
varying this value by a factor of 10 can produce only a factor of 
several variations in M™ n , owing to a rather week dependence 
of the total gas mass on «o,o ( see Fig. [3]). 

It is worth pointing out that the Kennicutt-Schmidt rela- 
tion ( TPTb is based on Ha measurements, which only reveal the 
presence of massive stars in the SF regions. Recently, mostly 
thanks to the GALEX satellite, measurements of UV fluxes be- 
came available for dwarf galaxies. It turns out that below ~ 1CT 2 
M yr _1 , the SFrate determined by the H a measurements largely 
underestimates that based on UV fluxes (Lee et al 1 12009b . If the 
IMF in DGs is top-light (i.e. steeper than the Salpeter slope and 
with a small upper m ass cutoff), Ho; fluxes can be ver y low even 
though SF is active dPflamm-Altenburg et al.L 12009b - This can 
affect the threshold value for star formation E c in equation ( fTTl i. 

9. Conclusions 

We have constructed a series of rotating equilibrium galaxies 
consisting of gas, stars, and a fixed DM halo, with masses of 
the latter M DM being in the 10 7 - 10 10 M range. Our mod- 
els differ from most previous studies in that we self-consistently 
take into account self-gravity of the gas component. Variations 
in the gas temperature, DM halo shape, rotation and non-thermal 
support, and also in the redshift for galaxy formation have been 
considered. We apply contemporary star formation criteria to the 
resulting equilibrium configurations to estimate the feasibility 
of large-scale star formation in our models. For the star forma- 
tion criteria, we choose the Toomre gravitational stability cri- 
terion with the Toomre parameter smaller than a critical value 
of Q c = 2.0 and the Kennicutt-Schmidt law with the gas sur- 
face density greater than a critical value of 2 cr i t = 5.0 M pc" 2 . 
These criteria need to be satisfied simultaneously at least in some 
parts of the gas disk in order for the model to be marked as star- 
formation-allowed (SFA). In addition, we require that the gas 
mass with number density greater than « = 1 .0 cm -3 exceed than 
10 4 M to allow for a SF event of non-negligible magnitude. 

We compare gas masses of the SFA models with the baryonic 
mass M b derived from the ACDM theory and WMAP4 data, for 
which the ratio of the baryon-to-DM mass is 0.17. We find the 
following: 

- For a given DM halo mass there exists a minimum gas mass 
M™ ln that is needed to achieve a state in which star forma- 
tion is allowed. The value of M™ m depends crucially on the 
gas temperature T g , the gas spin parameter a, the amount of 
nonthermal support in the gas disk, and, to a somewhat lesser 
extent, on the redshift for galaxy formation z s f- On the other 
hand, M™ n is rather insensitive to the form of the DM halo 
and to the pre-existing stellar disk, provided that the past SF 
efficiency does not exceed considerably 0. 1 . 

- As a rule, M™ n is smaller for galaxies with smaller a and T g 
and is greater for objects with greater non-thermal support 
o" e ff . In addition, M™ n may be smaller for objects that form 
at higher redshifts. 

- Depending on the gas temperature T g , gas spin parame- 
ter a, gas effective velocity dispersion cr eff , and the red- 
shift for galaxy formation z g t, the SFA models may have 
M™ n that is either greater or smaller than Mb. Models with 

M DM £ 10 9 M are usually characterized by M™ n 5 M b , 
implying that star formation in such galaxies is a natural 
outcome of their evolution. On the other hand, models with 
M DM ^ 10 9 M are often characterized by M™ n » M b , 



implying that they need much more gas than available ac- 
cording to the ACDM theory to achieve a state in which star 
formation is allowed. 

- A good agreement of our derived M™ n with Mb can be 
achieved if the spatially uniform gas temperature follows the 
virial relation T g oc M DM with a proper choice of the zero 
point or some added non-thermal support. However, the re- 
quired temperatures for objects with Mdm ^ 10 8 M are 
quite low (< a few hundred Kelvin) and the rotation curves 
are in poor agreement with those observed in DGs. 

- SFA models with M DM ^ 10 8 M and T g > a few x 10 4 K 
have M™ n that greatly exceed both Mb and Mdm, implying 
that some star-forming DGs may be baryon-dominated. 

We find that the gas volume density distribution of our model 
galaxies is crucially sensitive to the gas temperature, spin pa- 
rameter, and redshift of galaxy formation, implying a variety of 
possible equilibrium realizations for objects with the same DM 
halo mass. This means that the evolution of a dwarf galaxy may 
follow different pathways after the onset of star formation, de- 
pending on the values of T g , a, and z s f even for the same DM 
halo mass. 

Our modeling suggests that a star-formation-allowed state is 
more difficult to achieve in DM halos with mass ^ 10 9 M than 
in their upper-mass counterparts, because the required gas mass 
may be much greater than that available according to the ACDM 
theory. This implies that there may be a critical DM halo mass 
below which the likelihood of star formation and hence the total 
stellar mass drops substantially. Interestingly enoug h, the stel- 
lar ver sus DM halo mass relation recently derived by iGuo et all 
d2oToh using the SDSS survey and Millennium Simulations im- 
plies the existence of such a threshold value. On the other hand, 
DGs with the gas plus stellar mass greater than that of the DM 
halo are not unheard of and recent observations of the mass- 
to-ligh t ratios in Virgo Cluster dwarf ellipticals by jToloba et alj 
(1201 ll) and in gas-rich DGs by IS waters et al.l (1201 II) point to the 
existence of such objects. These observations, along with our re- 
sults, suggest that the ACDM paradigm is not universal and sig- 
nificant deviations from the corresponding Mb -Mdm trend are 
feasible. 

Our study outlines the importance of gas self-gravity (ne- 
glected in practically all hydrodynamical studies of isolated 
DGs) in building equilibrium galaxies. The main argument in 
favour of neglecting the gas self-gravity has been based on the 
assumption that the total gas mass is always much smaller than 
that of the DM halo. As we have demonstrated, this assump- 
tion may be grossly violated, particularly for low-mass DM ha- 
los with Mdm ^ 10 8 M . The reason is that the baryonic matter 
has to dominate the dark matter in objects with low DM halo 
masses in order to achieve the Kennicutt-Schmidt SF criterion. 
We emphasize that our results are strictly applicable to DM ha- 
los that have accumulated their mass reservoir quasi-statically 
and remain to be justified for object that have undergone a se- 
ries of violent mergers. At the same time, our main conclusions 
are not affected by moderate perturbations in a quasi-equilibrium 
state and reasonable variations in the adopted values of X cr i t and 

&. 
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tive and ghost grid zones with solid and dashed lines, respec- 
tively. The zone centers are marked with filled/open circles for 
active/ghost zones. The arrows indicate the backward finite- 
difference scheme used to discretize spatial derivatives. The in- 
nermost grid cell highlighted with a backslash palette refers to 
the seed value of the gas surface density «o,o- 
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Appendix A: Solving for the steady-state equations 



Figure I A. 1 1 shows the N x N computational mesh employed to 
discretize the steady-state equations ® and (0). The active zones 
are outlined with the solid lines, whereas the two rows of ghost 
zones (representing the reflecting boundary conditions along the 
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z- and r-axes) are marked with the dashed lines. The active/ghost 
zone centers are denoted with filled/open circles. 

A class of problems that does not require the knowledge of 
the gas density at the outer z and r boundaries (i.e., at N + 1 grid 
zones) can be solved using the following procedure. We use a 
first-order backward difference scheme (schematically shown by 
the arrows) to obtain a finite-difference representation of equa- 
tions (f3j) and (|4]i for the case with a spatially uniform cr g 
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where the indices i and j correspond to the z and r coordinate 
directions, respectively. We note that densities are defined at the 
zone centers while gravitational accelerations and velocities are 
defined at the corresponding zone interfaces. In order for this 
difference scheme to work, one needs to define the gas density 
at the ghost zones (which equal to those at the nearest active 
zones) and also the gas density at the innermost active zone (1,1) 
denoted in the paper as the seed density «o o. The corresponding 
zone is highlighted with the backslash palette in Figure [AT] 

With this choice of the discretization scheme and boundary 
conditions, one can notice that Pg and Pg are known for 

every value of pi'' 7 and the latter can be found by a fast forward 
substitution algorithm if one proceeds from left to right along 
the r-direction, starting from the bottom layer of zones and ad- 
vancing one horizontal layer after another in the direction of in- 
creasing z. 

Appendix B: Testing equilibrium configurations 

An important reliability check on the solution procedure is to 
test how our equilibrium configurations can be handled by time- 
dependent numerical hydrodynamics codes. If our steady-state 
models are correct, than a galaxy should stay in rotational equi- 
librium for at least 500 Myr, a typical time of interest when sim- 
ulating the effect of supernova explosions in DGs. 

To perform such a test, we use our time-dependent numerical 
hydrodynamics code employed earlier to study the effect of SN 
explosions in DGs in the local Universe and and large redshifts 
dVorobvov et al.L 12004 IVorobvov & Basil 120051: IVasiliev et al.L 
12008b . We intentionally turn off cooling and heating to avoid the 
system drifting out of equilibrium due to thermal effects. For the 
test, we use the reference model with Mum = 10 9 M e , a = 0.9, 
T g = 10 4 K, and «o,o = 5.0 cm" 3 . The gas surface density, Q 
parameter, and velocity profiles of this model are shown by the 
red solid lines in the second top row of Figure |2] 

Figure |BTJ shows the mean relative error (solid line) and the 
maximum relative error (dashed line) in the gas volume density 
p g (top) as a function of time t in our test model. The relative 
errors (in per cent) are calculated at every grid cell asPl 



Ap. 



IPg,/./0-Pg,v(0)l 
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(B.l) 



and demonstrate the degree to which our equilibrium is held by 
the code during the time evolution. The mean relative errors Ap s 



12 An extended definition of the relative error that takes into account 
a situation when the gas density declines with time, and is thus nor- 
malized to the current value of p g (t) rather than to the initial one p g (0), 
yields very similar results. 
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Fig. B.l. Maximum and mean relative errors (dashed and solid 
lines, respectively) in the gas volume density (top) as functions 
of time t for the reference test model described in Appendix B. 
The errors are calculated relative to the initial equilibrium con- 
figuration at t — Myr. The dash-dotted lines show the relative 
errors in the absence of gas self-gravity. 



are calculated by averaging the individual errors Ap gj! j over the 
entire computational grid. As one can see, the mean relative er- 
rors never exceed 0. 1 %, meaning the equilibrium is well pre- 
served globally. The maximum relative error never exceeds 3 % 
and is kept below 1% during most of the evolution. We note 
the maximum errors occur in dynamically unimportant regions 
near the axes at large m. The gas temperature shows essentially 
the same behaviour. This test convincingly proves the robustness 
and reliability of our solution procedure. 

To demonstrate the importance of gas self-gravity and to 
perform the final check on our self-gravitating equilibrium 
configurations, we artificially turn off the gas self-gravity in 
our time-dependent numerical hydrodynamics code. The dash- 
dotted lines in Figure IBTTI present the resulted mean relative er- 
ror. It is obvious that neglecting the gas self-gravity results in 
a complete destruction of the equilibrium state, with the mean 
relative errors exceeding 10 6 for the gas volume density! 
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